Fluid Elasticity Can Enable Propulsion at Low Reynolds Number 
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Conventionally, a microscopic particle that performs a reciprocal stroke cannot move 
through its environment. This is because at small scales, the response of simple 
Newtonian fluids is purely viscous and flows are time-reversible. We show that by 
contrast, fluid elasticity enables propulsion by reciprocal forcing that is otherwise 
impossible. We present experiments on rigid objects actuated reciprocally in viscous 
fluids, demonstrating for the first time a purely elastic propulsion set by the object's 
shape and boundary conditions. We describe two different artificial "swimmers" that 
experimentally realize this principle. 
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A striking feature of Newtonian viscous flows is that they can be time-reversible. 1 This 
feature is often referred to as kinematic reversibility and has important consequences, for 
example, in fluid transport in micro- and nano-fluidic devices,- mixing,- and self-propulsion 
of microorganisms.- For microorganisms swimming in simple liquids, linear viscous stresses 
that scale as pV/L are much larger than stresses from nonlinear fluid inertia, which scale 
as pV 2 , where V and L are the characteristic velocity and length-scale, and p and p are 
the fluid density and viscosity. For these swimmers, the ratio of inert ial to viscous stresses, 
calculated as the Reynolds number Re = pVL/p, is often 10 -3 or smaller. Because the 
swimmer has a density comparable to that of the fluid, its own inertia is also negligible. 
The resulting kinematic reversibility means that only non-reciprocal deformations of the 
swimmer can break time-reversal symmetry and result in net motion; this is known as the 
"scallop theorem." 5 

However, the hydrodynamic stresses on a microorganism need not be purely viscous. 
Many microorganisms live in complex fluid media that contain solids and/or polymers 6 . 
Fluids such as gels, mud, intestinal fluid, and human mucus are not Newtonian, and often 
possess viscoelastic behavior. Recent work has begun to explore the important higher-order 
effects of fluid elasticity on swimmers that can also move through Newtonian (non-elastic) 
fluid. Theoretical and numerical studies have shown that fluid elasticity can significantly af- 
fect the propulsion speed and efficiency of microorganisms,- - — and breaks the time-reversal 
symmetry between pushers and pullers. 15 Controlled experiments have shown that fluid elas- 
ticity usually hinders propulsion compared to Newtonian fluid, 16 although there is evidence 
of an increase in propulsion speed for rotating helices in highly elastic fluids.— It is becoming 
increasingly clear that the presence of elastic stresses in the medium can modify swimming 
in a nontrivial way. 

The possibility, however, that fluid elasticity can enable rather than modify propulsion, 
circumventing the scallop theorem, is still largely unexplored. Propulsion enabled by fluid 
elasticity has been predicted for three special cases of reciprocal motion: a flapping surface 
extending from a plane^ 1 ^; a sphere which generates small-amplitude sinusoidal motion of 
fluid along its surface 6 ; and a "wriggling" cylinder with reciprocal forward and backward 
strokes at different rates.— However, there remains no experimental demonstration, and such 
propulsion of free, finite-amplitude swimmers has not been studied at all. 

In this Letter, we consider the question of whether viscoelasticity alone can enable propul- 
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FIG. 1. (a) Typical "swimmers." Two epoxy beads are joined by steel wire to form polar (asym- 
metric) and symmetric dimers. (b) Top view of experiment. Two aligning electromagnets at 
constant current are orthogonal to two driving magnets, controlled by a computer, (c) The dimer 
with magnetization rh experiences torque r mag to align with the magnetic field, (d) Dimer geom- 
etry. Dimer orientation a oscillates around (a), which is parallel to 5 a iign- If a wall is present, the 
smallest separation between it and the dimer is the gap size h\ typically 30 fim at the start of the 
experiment. The length of the dimer is 2i?dimer; the bead at each end has radius i?bead- 

sion in the absence of inertia (i.e. low Re) by actuating a single rigid object reciprocally 
in a very viscous fluid. A propeller ("swimmer") such as one shown in Fig. [I] is immersed 
in a fluid and repeatedly reoriented by a magnetic field. The effects of inertia are absent 
due to high fluid viscosity (~4x 10 4 cSt), resulting in Re < 10 -4 , comparable to that of a 
swimming microorganism. By applying only magnetic torques, our apparatus reciprocally 
actuates just one degree of freedom in the system, the dimer's orientation a. For a purely 
viscous Newtonian fluid at low i?e, we find no net motion because a(t) is cyclic. Yet when 
a small amount of polymer is added to the fluid, making it viscoelastic, the same "stroke" 
results in propulsion, in a direction set by the dimer's shape and boundary conditions (cf. 
Fig. [2j). While the dimers are not strictly self-propelled and so are not true swimmers, the 
magnetic field provides only a reciprocal torque and does not itself create or direct propul- 
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FIG. 2. (color online) Locomotion at low Re due to fluid elasticity, (a) Polar dimer at De = 5.7. 
Left: dimer superimposed on silhouettes of motion over 1 cycle at /drive — 0.4 Hz. At top and 
bottom the dimer is placed in opposite orientations, but with the same actuation in each case. 
Middle, right: corresponding stroboscopic centroid trajectories, plotted over 50 s of driving at 
2.8 Hz in viscoelastic fluid, proceeding from the large symbol to the small one. In Newtonian fluid 
there is negligible net translation from e.g. inhomogeneities in the magnetic field. The direction in 
viscoelastic fluid is set by the dimer shape; in Newtonian fluid, it is not. (b) Net displacement vs. 
time for upper row of (a), with same symbols. The black "static" line shows motion with i? a ii gn 
only (no driving), (c), (d) Corresponding plots for a symmetric dimer next to a wall with 9 W = 45° 
and De = 0.8. As the dimer moves away from the wall, velocity perpendicular to the wall decays, 
while velocity parallel to the wall is constant. 

sion. This is thus the first experimental demonstration of purely elastic propulsion, wherein 
fluid elastic stresses are the sole source of net motion. 

Two experimental systems are used: (i) a polar (asymmetric) dimer far from any bound- 
aries and (ii) a symmetric dimer near a wall, as shown in Fig. [D(a). Each dimer consists of 
a piece of carbon steel wire of length 2i?dimer = 2.5-3 mm and diameter 230 /xm, with an 
epoxy bead of diameter 2i?bead ~ 500 fim at each end. The dimer has orientation a and 
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is magnetized with moment fa = am, so that a uniform magnetic field B reorients it with 
torque r mag = m x B, as depicted in Fig. [H(c). For experiments with the dimer next to a 
wall, a glass cover slip serves as a flat, vertical boundary (Fig. [T[cl]). 

The dimer is immersed in a container (50 mm tall, 30 mm in diameter) of either Newtonian 
or polymeric fluid (Fig.[T(b]). The Newtonian fluid is a 96%-corn syrup aqueous solution (by 
mass) with kinematic viscosity \ij p of approximately 4 x 10 4 cSt. The polymeric solution 
is made by adding 0.17% (by mass) of high-molecular- weight polyacrylamide (PAA, M w = 
10 6 ) to a viscous Newtonian solvent (93%-corn syrup aqueous solution). The solution is 
considered dilute: the overlap concentration c* for PAA is ~ 0.34% (c/c* = 0.5).— Using 
a strain-controlled rheometer, we find that the PAA solution is an elastic fluid with nearly 
constant viscosity, varying with a power law index of n — 0.96 up to 50 s -1 . 21 This way, the 
effects of shear-thinning viscosity can be decoupled from those of elasticity.—^ The fluid 
relaxation time A, measured in stress relaxation tests, is approximately 2 s. 21 

A schematic of the apparatus is shown in Fig. [H(b). Four electromagnets reorient the 
dimer in a fluid cell. A reciprocal "wiggling" motion is achieved with two diametrically 
opposed electromagnets generating a constant field i? a i ign , and a pair of orthogonal magnets 
generating the AC field B& r - lYe . The driving magnet current is controlled by a computer via 
a power amplifier. The magnitude of £? a i ign and amplitude of B& r - lYe are (9(10 3 G). The 
amplitude of dimer rotation is nearly 45°, decaying at high /drive in some experiments as 
detailed below. For polymeric fluid, we define the Deborah number De, the product of the 
longest fluid relaxation time A and the driving frequency /drive- 
in addition to the magnetic torque, the dimer also experiences an undesirable translational 
force 

F mag = V(m • B) (1) 

due to inhomogeneities in the magnetic field. The dimer is positioned so that when both 
the aligning and driving fields are at full strength, the translation velocity from F mag is 
< 0.3 /im/s. This velocity is an upper bound on the translation in our experiments that can 
be attributed to F mag : the sinusoidally- varied driving field is only briefly at full strength, 
and fa and B are not aligned as the dimer rotates, reducing the dot product in Eq. ED 

Evidence of purely elastic propulsion is shown in Figs. [2](a, b) for the polar dimer far 
away from boundaries at De = 5.7 (/drive = 2.8 Hz; Re = 1.2 x 10 -4 ) and Fig. El(c,d) for the 
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symmetric dimer near a wall at De = 0.80 (/drive — 0.4 Hz; Re = 6.6 x 10 -5 ). The dimers 
are imaged using a CCD camera to extract orientation and centroid position. The camera 
is aligned with the vertical axis and the apparent horizontal motion from sedimentation is 
< 0.05 /im/s. The data show a striking contrast between performing reciprocal motion in 
Newtonian and in viscoelastic fluid. Figure [2](a, b) shows that in viscoelastic fluid, far from 
any boundaries, the polar dimer is able to achieve net motion at constant speed even under 
reciprocal forcing; no net motion is observed in a Newtonian fluid under the same conditions. 
The polar dimer moves in the direction of its large end, as shown by the arrows in Fig. [2](a). 

Net motion is also observed for a symmetric dimer near a wall in viscoelastic fluid 
(Figs. [2[c, d]). The symmetric dimer translates away from and along the wall with ap- 
proximate direction (a), as shown in Fig.[2](c). Corresponding behaviors are seen in variants 
of the geometry where the dimer is flipped by 180° or 9 W is varied in increments of 90°. 
Figure [2](d) shows the typical net displacement of these trajectories as a function of time. 
After a short transient, the symmetric dimer immersed in a viscoelastic fluid achieves a con- 
stant velocity, primarily parallel to the wall. This case is representative of the behavior for 
a wide range of Q W1 excluding the limiting cases of 9 W ~ 90° (near-negligible propulsion) and 
@w ~ 0° (motion primarily away from wall). This symmetry-breaking is distinct from that 
experienced by a conventional low-i?e swimmer in Newtonian fluid, wherein a wall alters 
the trajectory of a swimmer or particle, but does not change the fundamental nature of 
propulsion^i^; there, the particle can self-propel without the wall. 

It important to note that in a Newtonian fluid, all experiments discussed above yield 
negligible net displacement, comparable to the effects of F mag and sedimentation when driv- 
ing is turned off altogether (Fig. [2[b, d]). Furthermore, the direction of displacement in 
Newtonian fluid is not controlled by particle shape or boundary conditions, confirming that 
it is not hydrodynamic in origin. Experiments with a symmetric dimer in a viscoelastic 
medium far from any boundaries also yield negligible net displacement (not shown). 

Our experiments show that the net motion achieved by the dimers in the polymeric solu- 
tion results from elastic stresses due to flow-induced changes in polymer conformation. These 
elastic stresses are history-dependent and do not entirely cancel out over one forcing period, 
but instead have a small rectified component that accumulates. A rheological property of 
polymeric solutions that is of particular relevance here is the first normal stress difference 
Ni = tqo — r rr , where r, 9 and z are cylindrical coordinates and r is the fluid stress tensor. 
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FIG. 3. (color online) Instantaneous streamlines, computed from particle tracking experiments, 
as a polar dimer rotates far from any boundaries, here plotted in the dimer frame. The different 
streamline curvature at each end suggests that forces due to the fluid's normal stress difference Ni 
are unbalanced, contributing to propulsion. 

Ni grows nonlinear ly with fluid strain rate and, to lowest order, scales with strain rate as 
7 2 ,— consistent with measurements of our own polymeric fluid.— The combination of Ni and 
curved streamlines in a given flow results in an inward-pointing volume force —Ni/r in the 
radial direction. In Fig. [3] we show instantaneous streamlines during the stroke of the po- 
lar dimer, computed from experimentally measured velocity fields. The curved streamlines 
around each bead, and the asymmetry in that curvature due to the dimer shape, strongly 
suggest that Ni plays a role in propulsion. The greater strain rate and curvature at the 
small end of the dimer suggest that a stronger volume force there will move the dimer in 
the direction of its large end, consistent with observed propulsion. 

To gain further insight into possible mechanisms, the effects of elasticity on propulsive 
speed are investigated for the polar and symmetric dimers, as shown in Fig. [H The impor- 
tance of elasticity is quantified by De. For both geometries, translation speed increases with 
elastic stresses for the range of De investigated here. This trend is also seen for comparable 
De in the cases of purely elastic propulsion numerically and theoretically investigated to 
dateiUU&U 

As described above, the different streamline curvature at each end of the polar dimer sug- 
gests a role for normal stress difference effects as well as possible hydrodynamic interactions 
between the ends of the dimer. However, a description of elastic effects based on steady 
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FIG. 4. (color online) Dependence of mean propulsion on driving parameters. Both types of 
particles show an increase in propulsion with Deborah number De, consistent with an elastic 
phenomenon, (a) Polar dimer. As /drive is varied, propulsion appears to be controlled by De, 
not Wi (inset), suggesting that the mechanism is not suited to a steady- flow description. Here 
7 = ^i?dimer/^bead 5 where 9 is angular velocity. Note that strain amplitude decays as /drive is 
increased, due to viscous resistance, accounting for the difference between the figure and its inset. 
This decay is the likely cause of the turnover at the highest /drive- (b) Polar dimer in a separate 
experiment (fluid viscosity ~ 40% higher) showing propulsion as a function of strain amplitude 
and Wi (inset). Magnetic current amplitude is varied at De = 2 (/drive — 0.7 Hz, red squares) 
and De = 4 (1.4 Hz, blue circles). Data show a much weaker dependence on strain amplitude 
and Wi than on De. (c) Symmetric dimer at wall, showing much greater velocities and different 
scaling. Velocity component parallel to the wall is plotted. The range of /drive is limited in order 
to maintain constant strain amplitude. Velocity measurements are cut off by noise at low De. 

shear, and thus on a single strain rate, is likely inadequate to model propulsion. This is seen 
in the insets of Fig. [4](a, b), where we plot velocity vs. Weissenberg number Wi = A7, using 
the characteristic strain rate \J (7 2 ). In experiments, we find that \J ( A / 2 ) is proportional 
to /drive times strain amplitude. Propulsion is much more sensitive to /drive than to strain 
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amplitude, and so is poorly characterized by Wi. The dynamics would therefore best be 
modeled by fully accounting for the unsteadiness of the flow. 

For the symmetric dimer results in Fig.[H(c), following the observations in Fig.[2](c, d), we 
choose to plot the dominant velocity component, parallel to the wall, that is approximately 
constant over each movie. To keep geometry constant with respect to the fixed wall, we keep 
strain amplitude nearly constant by limiting /drive- In this geometry, elastic effects increase 
nearly linearly with De, unlike the roughly quadratic case of the polar dimer. We also find 
that even at these lower /drive, translation velocity is an order of magnitude greater than 
that of the polar case. 

Propulsion of the symmetric dimer is inconsistent with two models we discuss here. First, 
while it is known that a particle moving steadily near a wall in viscoelastic fluid will experi- 
ence a lift force away from the wall that scales as Rbe&d/h (e.g.—), Figs.[2](c, d) show that the 
dominant motion we observe is parallel to the wall and decays little as h increases. Second, 
Fu et al. 12 show that for one type of swimmer in viscoelastic fluid, performing different parts 
of a reciprocal stroke at different rates enables propulsion. In our experiment, the dimer 
tip approaches the wall < 10% slower than on the return stroke, due to the nature of the 
magnetic driving. But on this basis, the cited analysis suggests net motion in the opposite 
direction of what we observe. The elastic mechanisms for this translation at a wall, and the 
propulsion of the polar dimer, thus remain open questions. 

In summary, we have made the first experimental demonstration of propulsion by fluid 
elasticity, using reciprocally actuated dimers at low Re. The actuation yields no net motion 
in Newtonian fluid, but a small amount of polymer in the fluid adds an elastic response 
to the driving, permitting propulsion. The dimer shape directs propulsion in elastic fluid, 
without nearby boundaries. Near a wall, a dimer is propelled both parallel to and away 
from the boundary. All other propulsive strategies at low Re in linear (e.g. Newtonian) 
fluids require a non-reciprocal stroke. Here, a swimmer may employ a reciprocal stroke 
provided the fluid is viscoelastic. We note that because time-reversibility is broken by the 
fluid and not the stroke, a time-reversal of driving does not reverse propulsion, as it would 
for a non-reciprocal swimmer; here, the direction of propulsion is set by geometry alone. 

We estimate propulsive (Froude) efficiency, defined as the fraction of total hydrodynamic 
power corresponding to net motion, to be (9(1%) — comparable to the measured ~ 2% effi- 
ciency of the common low-i?e swimmer E. Coli— The relatively slow propulsion we observe, 
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at most ~ 5 /xm/s, is due to the high fluid viscosity required to make sedimentation and 
shear-thinning effects negligible. We expect significantly faster propulsion in other realiza- 
tions of this principle particularly for smaller dimers which would allow for viscosities that are 
much smaller. Also, fluid elasticity effects are expected to become more pronounced as the 
dimer (or swimmer) is miniaturized since the elasticity number, defined as El = A/i/ pL 2 ) 
scales inversely with the square of the dimer length scale L. Finally, at higher De or in 
geometries with greater streamline curvature, purely elastic instabilities may cause sponta- 
neous propulsion without a wall or an asymmetric dimer shape — or may greatly enhance 
the mechanisms demonstrated here.—"— While organisms may not exploit the principle de- 
scribed here, viscoelastic media are common in nature, 6 and reliance on nonlinear rheology 
is not without precedent.—^ 

Our work is also a proof-of-concept for an artificial "swimmer" that moves through com- 
plex fluid with only reciprocal actuation, a simple body shape, and no moving parts — a 
less complicated design than for other propulsive strategies.—^ These principles could also 
be applied to pumps,— ^ or to exploiting other types of nonlinear fluid rheology. Further 
understanding of this effect and similar ones could greatly simplify fabrication of micro- 
swimmers for many artificial environments, or for biological settings where viscoelasticity is 
ubiquitous. 
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